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Abstract — Numerical  simulations  are  used  to  solve  the  coupled 
partial  differential  equations  describing  stimulated  Brillouin 
scattering  (SBS)  built  up  from  random  thermal  phonons  as  a 
function  of  time  and  the  longitudinal  spatial  coordinate  in  an 
optical  fiber.  In  the  case  of  a  passive  fiber,  a  laser  beam  is  incident 
with  constant  power,  but  its  frequency  is  linearly  ramped  at 
1.55  fL m  at  a  rate  of  up  to  1016  Hz/s.  High  chirp  rates  lead  to 
an  increased  Brillouin  spectral  bandwidth  and  decreased  gain. 
The  resulting  SBS  suppression  is  well  described  by  an  adiabatic 
model  and  agrees  with  experimental  results.  For  an  18-m  active 
fiber  pumped  at  1.06  p m  and  chirped  at  up  to  2  x  1016  Hz/s, 
the  suppression  enables  output  laser  powers  in  the  kilowatt  range 
while  maintaining  a  narrow  instantaneous  linewidth. 

Index  Terms — Brillouin  scattering,  chirped  lasers,  fiber 
amplifiers,  numerical  simulation 


I.  Introduction 

BRILLOUIN  scattering  is  one  of  the  lowest  order  nonlinear 
effects  that  arises  in  optical  fibers  and  it  thus  limits  the 
transmitted  laser  power  P^.  Above  a  threshold  incident  laser 
power  Pth,  the  Brillouin  scattering  becomes  stimulated  rather 
than  spontaneous  and  the  Stokes  backscattered  power  Ps  rises 
dramatically.  Various  methods  have  been  proposed  to  increase 
the  threshold  [1]— [5].  In  the  present  paper,  we  analyze  the 
effect  of  linearly  sweeping  the  frequency  of  the  pump  laser  at 
a  fast  chirp  rate  /?.  To  be  effective,  the  chirp  must  be  large 
enough  that  the  pump  laser  is  swept  out  of  the  Brillouin  gain 
bandwidth  Avb  within  the  transit  time  of  the  fiber,  tf  =  nL/c, 
where  L  is  the  length  of  the  fiber,  n  is  the  core  refractive  index, 
and  c  is  the  speed  of  light.  That  is,  one  requires  /I  to  be  large 
enough  that  tf  >>  tc  where  tc  =  Avb/2/?.  (The  factor  of 
2  comes  from  the  fact  that  the  laser  and  Stokes  waves  are 
counterpropagating  in  the  fiber.)  Here  we  analyze  values  of  /? 
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up  to  2  x  10 16  Hz/s,  close  to  the  upper  range  of  chirping  that 
has  been  demonstrated  to  date  for  a  diode  laser  by  ramping 
its  current  [6],  [7]. 

Section  II  describes  numerical  simulations  of  the  Brillouin 
scattering  in  a  passive  fiber  with  an  incident  laser  beam  that 
is  chirped.  Section  III  then  presents  the  simulation  results 
for  a  6-km  fiber.  Such  long  fibers  have  a  low  threshold 
because  Pth  ~  2\A/goL  in  the  unchirped  case,  where  A  is 
the  modal  area  and  go  is  the  peak  Brillouin  gain  coefficient. 
Consequently,  low  laser  powers  can  be  used  to  experimentally 
verify  the  simulations.  Agreement  is  found  between  theory 
and  measurements.  Next,  Sec  IV  presents  results  of  the  time- 
dependent  simulations  for  a  short  fiber  (L  =  17.5  m)  having 
characteristics  similar  to  those  used  for  high-power  laser 
delivery.  Correspondingly  larger  chirps  are  needed  to  suppress 
stimulated  Brillouin  scattering  (SBS)  in  this  case.  The  results 
agree  with  an  adiabatic  model  that  depends  on  frequency 
(instead  of  time)  for  chirps  ranging  up  to  10 16  Hz/s.  Finally, 
Sec  V  adds  laser  gain  to  the  simulations  to  predict  the  output 
power  from  a  fiber  amplifier.  Chirping  the  input  laser  beam 
at  2  x  10 16  Hz/s  enables  one  to  reach  the  kilowatt  range  by 
increasing  the  SBS  threshold. 


II.  Theory  for  the  Dynamic  Simulations 


The  Brillouin  scattering  is  modeled  by  three  complex  cou¬ 
pled  partial  differential  equations  (PDEs)  [8], 


d£L  n  d£L  C  ~ a 
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We  checked  that  adding  the  two  extra  derivative  terms  from  (8) 
in  [5]  to  the  left-hand  side  of  (lc)  did  not  noticeably  change 
our  simulation  results;  so  for  computational  simplicity,  those 
terms  are  omitted.  Here  the  laser  electric  field  £l(z, 
Brillouin  Stokes- shifted  field  Es(z,t),  and  density  variation 
p(z,t)  of  the  fiber  from  its  mean  value  po  depend  on  time  t 
and  longitudinal  position  z,  varying  from  z  =  0  at  the  front 
face  of  the  fiber  to  z  =  L  at  the  rear  face.  The  spectral  full- 
width-at-half-maximum  (FWHM)  of  the  spontaneous  Brillouin 
peak  is  Avb.  For  a  silica  fiber  at  an  incident  laser  wavelength 
of  2l  =  1-55  pm,  the  refractive  index  is  n  =  1.447,  the 
mean  density  is  po  =  2210  kg/m3,  and  the  speed  of  sound 
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is  v  =  5960  m/s  [8].  The  loss  coefficient  in  the  fiber  is 
a  =  0.2  dB/km  =  0.0461  /km.  The  laser  gain  £*  is  taken  to 
be  independent  of  t  and  z  for  simplicity;  £*  =  0  for  a  passive 
fiber.  The  optic  coupling  parameter  is  [8] 


n  y 

2npol]^M 


(2) 


where  the  electrostriction  coefficient  is  y  —  0.902  [9]  and  the 
polarization  is  presumed  to  be  completely  scrambled  in  the 
fiber  so  that  M  =  1.5  [10].  The  acoustic  coupling  parameter 
is  [8] 


nnsoy 


(3) 


where  £o  is  the  permittivity  of  free  space.  The  Langevin 
noise  source  f(z,t)  is  delta-correlated  in  time  and  space  such 
that  [11] 


(f(z,  t)f  *  (z',  O)  =  Q  S(z  ~  z')d(t  -  O  (4) 


where  the  thermal  phonons  are  described  by  the  strength 
parameter 

AnkTpo  Avb 

2  =  — (5) 

Here  k  is  Boltzmann’s  constant,  T  =  293  K  is  room  temper¬ 
ature,  and  A  is  the  fiber  modal  area. 

The  three  PDEs  (1)  are  solved  by  iterated  finite-difference 
approximations  on  a  grid  of  time  and  space  points  [12].  The 
spatial  step  size  dz  is  chosen  to  be  small  enough  that  a  factor 
of  two  reduction  does  not  change  the  final  values  in  (8)  below 
by  more  than  10%.  The  temporal  step  size  is  dt  =  ndz/c. 
The  boundary  conditions  are  Es(L,t )  =0  and 


£l(0,0 


(6) 


where  Po  is  the  constant  laser  power  incident  at  z  =  0.  The 
argument  of  the  exponential  is  the  chirped  phase 


(pit,  z)  =  2n 


t+n(L-z)/c 

J  fit'dt' 

o 


(7) 


evaluated  at  z  =  0.  The  three  complex  fields  are  first  calculated 
along  the  fiber  for  times  up  to  t  =  20/f  to  ensure  that  the 
initial  relaxation  oscillations  [13]  have  decayed  away.  The 
equations  are  then  iterated  over  5  more  transit  times  to  acquire 
a  statistical  average.  The  transmitted  laser  power  P and 
reflected  Stokes  power  Ps  are  computed  by  averaging  over 
those  additional  steps 


PL  =  2nc£QA[\E^iL,  t)\2)  and 

Ps  =  jnce0A(\Es(0,t)\2).  (8) 

The  results  were  checked  by  verifying  that  Pl  ~  (Po  —  Ps) 
exp  (— aL )  in  a  passive  fiber,  and  agreement  was  found  to 
within  about  1%  at  zero  chirp. 


Fig.  1.  Brillouin  reflectivity  versus  power  for  both  an  unchirped  (/?  =  0) 
and  a  chirped  (at  ft  =  10 14  Hz/s)  laser  beam  incident  on  a  long  fiber. 


III.  Dynamic  Results  for  a  Long  Passive 
Fiber  Compared  to  Experiment 

We  calculated  Ps  as  a  function  of  Po  f or  a  passive  6-km 
single-mode  fiber.  The  FWHM  of  the  spontaneous  Brillouin 
peak  was  measured  to  be  Avb  =39  MHz  using  a  heterodyne 
technique  with  an  RF  spectrum  analyzer  [7].  The  fiber  modal 
area  was  measured  to  be  A  =  nr2  with  r  =  4.55  pm  by  a 
knife-edge  scan  in  the  far  field. 

The  reflectivity  Ps/Po  is  plotted  in  Fig.  1  for  both  an 
unchirped  and  a  chirped  laser.  For  incident  powers  well  below 
threshold,  the  reflectivity  levels  off  to  a  spontaneous  value 
that  is  independent  of  chirp  and  is  in  good  agreement  with 
the  experimental  value  of  Po  =  (3.0  zb  0.5)  x  10-6  measured 
using  an  optical  spectrum  analyzer  (OS A).  On  the  other  hand, 
for  large  incident  powers  the  reflectivity  approaches  100%. 
The  threshold  (defined  as  the  incident  power  at  which  the 
reflectivity  is  equal  to  1%)  is  approximately  two  orders  of 
magnitude  larger  in  the  chirped  case  than  in  the  unchirped 
case. 

To  directly  compare  these  results  to  the  experimental  data, 
the  total  backscattered  power  is  next  plotted  in  Fig.  2  as  a 
function  of  Po-  This  backscattered  power  is  the  sum  of  the 
Rayleigh  power  Pr  =  2.3  x  10_4Po  determined  from  low- 
power  OSA  measurements  and  the  Brillouin  Stokes  power  Ps 
from  Fig.  1.  The  agreement  is  reasonable  between  experiment 
and  theory  with  no  free  parameters. 

The  peak  value  of  the  Brillouin  gain  coefficient  is 


6U  -2  a  A/r  V  ' 

ncpov  a^AvbM 

which  equals  6.4  x  10“ 12  m/W  for  the  parameters  given 
above.  This  value  is  a  factor  of  2  to  4  times  smaller  than 
what  is  typically  measured  for  fibers  [14],  possibly  because  of 
inhomogeneous  broadening  of  the  Brillouin  peak.  According 
to  (9),  it  is  only  the  product  £oAvb  that  should  be  the  same 
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incident  laser  power  (W) 


chirp  (Hz/s) 


Fig.  2.  Comparison  of  the  simulated  and  experimental  [7]  total  backscattered 
power  for  the  unchirped  and  chirped  laser  sources. 


for  all  polarization- scrambling  silica  fibers  pumped  at  1.55 
and  not  their  individual  two  values. 

The  power  spectral  density  (in  dBm/Hz)  of  the  reflected 
Stokes  electric  field  Es(0,t)  was  plotted  for  each  simulation 
run  as  a  function  of  frequency  up  to  the  Nyquist  limit.  Line- 
shape  functions  were  then  fit  to  these  spectra.  The  resulting 
FWHM  are  plotted  as  the  dots  in  Fig.  3  for  the  long  fiber  at 
Po  =  1  mW,  well  below  threshold,  for  four  different  chirps.  At 
the  two  lower  chirps  of  10 10  and  10 12  Hz/s  the  lineshapes 
were  found  to  be  Lorentzian,  whereas  at  the  two  higher  chirps 
of  5  x  10 12  and  10 13  Hz/s  the  lineshapes  were  Gaussian.  The 
continuous  curve  is  a  plot  of  the  broadening  expected  from  a 
simple  model,  namely 

Ay  =  Avb  +nLfi/c  (10) 

which  fits  the  simulation  results  to  within  the  error  bars  on 
the  widths.  Replacing  Aub  by  Ay  in  (9),  this  increase  in  the 
Brillouin  bandwidth  implies  a  reduction  in  the  peak  gain  and 
hence  a  higher  threshold  Pth  compared  to  the  unchirped  case. 

IV.  Dynamic  Results  for  a  Short  Passive 
Fiber  Compared  to  Other  Models 


Fig.  3.  FWHM  Av  of  the  Brillouin  line  as  a  function  of  /?. 


Fig.  4.  Brillouin  reflectivity  for  a  laser  beam  with  a  linear  chirp  /?  ranging 
between  10 10  and  1016  Hz/s  that  is  incident  on  a  short  fiber. 


should  be 


2  n2kTy2La 
Mn2  pov2  A^A 
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Consider  a  passive  delivery  fiber  having  a  length  of  L  = 
17.5  m,  a  core  radius  of  r  =  13.75  jum ,  and  a  spontaneous 
Brillouin  FWHM  of  Apb  =  20  MHz  conservatively  chosen 
to  be  on  the  low  end  of  the  range  of  linewidths  measured  for 
silica  fibers  [15]  so  that  the  peak  Brillouin  gain  is  as  large 
as  it  might  be  (namely  go  =  1.2  x  10-11  m/W).  All  other 
parameters  are  taken  to  be  the  same  as  those  listed  in  Secs  II 
and  III.  The  reflectivity  is  plotted  in  Fig.  4  for  seven  different 
values  of  the  chirp  /?. 

Well  below  threshold,  the  reflectivity  has  a  constant  value 
of  Ro  ~  5  x  10“ 10  independent  of  chirp  and  pump  power. 
A  heuristic  steady- state  model  [8]  for  the  buildup  of  the 
spontaneous  Stokes  wave  from  thermal  noise  predicts  that  it 


where  the  effective  absorption  length  of  the  fiber  is  La  = 
[1  —  exp(— aL)\/a.  This  model  implies  Ro  ~  3  x  10  9  which 
is  only  a  factor  of  6  larger  than  our  simulation  value.  Similar 
agreement  is  found  for  the  long  fiber  in  Sec  III:  (11)  predicts 
~  9  x  10  6  whereas  Fig.  1  has  a  low-power  reflectivity  of 
i?0~2x  10“6. 

The  curves  in  Fig.  4  cross  1%  reflectivity  at  the  threshold 
incident  laser  powers,  /V  The  interpolated  crossing  points 
are  plotted  in  Fig.  5  as  the  squares.  At  the  maximum  chirp 
of  10 16  Hz/s,  the  threshold  has  increased  by  a  factor  of  50 
compared  to  an  unchirped  laser  source. 

In  Fig.  5,  the  threshold  power  scales  linearly  with  the 
chirp  /?  above  10 14  Hz/s.  These  results  agree  with  an  adiabatic 
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chirp  (Hz/s) 


Fig.  5.  Threshold  power  Pt^  as  a  function  of  the  chirp  /?.  The  squares  are 
from  the  time-dependent  simulations  and  the  continuous  curve  is  based  on  an 
adiabatic  model  [7]. 


Fig.  6.  Transmitted  laser  power  as  a  function  of  incident  laser  power  for 
various  chirps. 


model,  according  to  which  the  threshold  should  increase  as 
tf/tc  where  the  fiber  transit  time  is  tf  =  nL/c  and  the  time 
required  for  the  pump  laser  to  chirp  out  of  the  Stokes  band¬ 
width  is  tc  =  Avb/2/?.  Combining  a  resonant  integration  [7] 
of  this  increase  with  the  familiar  [16]  factor  of  21  for  the 
unchirped  threshold  value  goP±L/A,  one  obtains 


21A  tf/tc 

Pth  = - T— -  (12) 

goL  tan  '(tf/tc) 

plotted  as  the  continuous  line  in  Fig.  5.  The  agreement  with 
the  time-domain  simulations  is  remarkable.  Furthermore,  one 
can  verify  graphically  that 


tf/t  C  ^  1  tf  1 

tan_1(ff/  tc)  2  tc  Avb 


Avb  + 


nLp- 

c 


(13) 


thereby  showing  that  (10)  is  also  consistent  with  the  adiabatic 
model. 


V.  Predicted  Output  Laser  Power 
for  a  Short  Active  Fiber 

Chirping  enables  one  to  maintain  gain  in  a  fiber  amplifier 
by  suppressing  the  power  lost  to  Brillouin  backscattering.  To 
demonstrate  this  effect,  we  simulated  an  ytterbium-cladding- 
doped  fiber  with  length  L  =  18  m  and  modal  radius 

r  =  13.75  /im  at  a  laser  wavelength  of  2l  =  1.064  jum 
with  a  spontaneous  Brillouin  FWHM  of  Avb  =  50  MHz 
representing  a  round  average  of  typical  linewidths  reported 
in  the  literature  [16].  As  a  simple,  conservative  model  of 
10  x  amplification,  the  gain  coefficient  for  both  the  laser  and 
Stokes-shifted  Brillouin  waves  was  set  equal  to  £*  =  ln(10)/L 
in  (1).  (A  more  accurate  model  would  simultaneously  solve 
the  rate  equations  for  the  populations  of  the  ground  and 
excited  laser  levels  to  compute  the  gain  from  a  separate  pump 
source  [17].)  In  Fig.  6,  values  of  Pl  calculated  from  (8) 
are  plotted  versus  Po  for  chirps  /?  between  2  x  1014  and 
2  x  10 16  Hz/s.  At  low  chirps,  the  transmitted  power  rolls 


off  when  it  surpasses  threshold.  For  example,  at  a  chirp  of 
2  x  1014  Hz/s  and  incident  laser  power  of  100  W,  which  is 
well  beyond  threshold,  in  the  absence  of  gain  the  laser  power 
drops  to  and  levels  off  at  87  W  within  the  first  few  meters  of 
the  fiber,  while  the  Stokes  power  correspondingly  rises  from 
0  to  13  W  over  the  same  distance.  In  contrast,  with  the  10  x 
gain  switched  on,  the  laser  power  drops  to  a  minimum  after 
about  1  m  of  transit  but  is  then  amplified  over  the  remainder 
of  the  fiber’s  length  to  215  W,  as  plotted  in  Fig.  6;  meanwhile 
the  Stokes  beam  begins  to  grow  a  few  meters  from  the  rear 
end  of  the  fiber  and  attains  about  400  W  at  the  front  face.  At 
the  highest  chirp,  in  contrast,  10  x  amplification  is  maintained 
for  incident  powers  up  to  200  W  in  Fig.  6,  resulting  in  kilowatt 
output  laser  powers. 


VI.  Conclusion 

Time-dependent  numerical  simulations  show  that  SBS  in 
optical  fibers  can  be  suppressed  by  linearly  chirping  the  laser 
frequency.  One  significant  advantage  of  this  method  compared 
to  competing  techniques  that  broaden  the  incident  laser  beam 
with  white  noise  is  that  a  narrow  instantaneous  linewidth  is 
maintained,  enabling  coherent  combination  of  the  outputs  of 
several  fiber  lasers  for  power  scaling.  Linear  chirping  allows  a 
mismatch  in  the  length  of  one  amplifier  relative  to  another  to 
be  compensated  by  acousto-optically  shifting  the  frequencies 
at  the  input  [18]. 

These  dynamic  simulations  are  in  good  agreement  with 
experimental  results  for  a  6-km  fiber.  The  simulations  should 
therefore  also  be  valid  for  a  short  17.5-m  fiber  having  char¬ 
acteristics  appropriate  to  high-power  laser  systems.  The  time- 
dependent  numerical  results  for  a  short  passive  fiber  are  in 
good  agreement  with  an  analytic  adiabatic  model,  thereby 
validating  use  of  that  model  to  predict  the  performance  of 
chirped  fiber  lasers  without  requiring  the  demanding  computer 
overhead  of  the  full  time-dependent  calculations.  The  results 
for  a  short  active  fiber  indicate  that  laboratory  attainable 
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chirps  enable  amplification  to  kilowatt  output  powers  while 
maintaining  a  narrow  instantaneous  laser  linewidth. 
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